Defect-mediated turbulence in systems with local deterministic chaos 
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We show that defect-mediated turbulence can exist in media where the underlying local dynamics 
is deterministically chaotic. While many of the characteristics of defect-mediated turbulence, such 
as the exponential decay of correlations and a squared Poissonian distribution for the number of 
defects, are identical to those seen in oscillatory media, the fluctuations in the number of defects 
differ significantly. The power spectra suggest the existence of underlying correlations that lead to 
a different and non-universal scaling structure in chaotic media. 
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Weakly driven, dissipative, pattern- forming systems 
often exhibit the spatiotemporal chaos in the form of 
defect-mediated turbulence, where the dynamics of a pat- 
tern is dominated by the perpetual motion, nucleation 
and annihilation of point defects (or dislocations) 1]. Ex- 
amples can be found in striped patterns in wind driven 
sand, electroconvection in liquid crystals |2|, nonlinear 
optics , fluid convection 0, Q , autocatalytic chemical 
reactions and Langmuir circulation in the oceans Q. 
These results suggest that the dynamics of these very 
different systems can be characterized by a universal de- 
scription which is based only on the defect dynamics. 

A statistical description of defect mediated turbulence 
based on a simple model for the defect dynamics was 
given by Gil, et al. 8]. Treating the defect pairs as statis- 
tically independent entities, the nucleation rate for pairs 
of defects was taken to be independent of the number of 
pairs n and, based on the topological nature of the de- 
fects, the annihilation rate was taken to be proportional 
to n 2 . This directly led to a squared Poissonian distri- 
bution for the probability distribution function (PDF) of 
n. Gil, et al. found that simulations of a spatiotempo- 
ral chaotic state of the complex Ginzburg-Landau equa- 
tion (CGLE), which is the prototype of oscillatory me- 
dia, agreed with their prediction. Rehberg, et al. Q 
measured the PDF of n for defect-mediated turbulence 
in electroconvection of nematic liquid crystals and found 
it to be consistent with the predicted squared Poissonian 
distribution. Later, Ramazza, et al. 3] investigated a 
defect turbulent state in optical patterns and found that 
their data were not conclusive. Very recently, Daniels 
and Bodenschatz 5] studied the defect-mediated turbu- 
lent state of undulation chaos in inclined layer convection 
of a fluid and found that the observed pair nucleation and 
annihilation rates agree with the theoretical predictions 
in 0. They also derived the PDF of n for the case that 
single defects can leave and enter through the boundaries 
and found agreement with their measurements. When 
boundary effects are negligible, this PDF reduces to the 
squared Poissonian distribution. 

The experimental and theoretical studies of defect- 



mediated turbulence have focused exclusively on first or- 
der statistics like the PDF of n and creation and annihila- 
tion rates. Moreover, theoretical investigations have been 
carried out only for media with underlying oscillatory dy- 
namics, generically described by the CGLE. Yet, in many 
cases the local dynamics may differ from simple oscilla- 
tory behavior; instead, complex-periodic or even chaotic 
attractors may exist (see, e.g., Ref. 0). Chemically re- 
acting systems, notably the Belousov-Zhabotinsky (BZ) 
reaction, are known to exhibit deterministic chaos ml. 
This leads to a variety of new spatiotemporal states [ll| . 
In this Letter, we show that defect-mediated turbulence 
can exist in media where the underlying local dynamics 
is chaotic and that second-order statistics can be used to 
distinguish between different media. This implies that a 
universal description of defect-mediated turbulence can- 
not encompass second-order or higher correlations. 

Our focus is on systems where the dynamics of the 
spatially homogeneous system, described by ordinary dif- 
ferential equations, has a deterministic chaotic attr actor; 
hence, at least three phase space variables are required 
in contrast to the two- variable descriptions of simple os- 
cillatory media. A specific example of such a chaotic sys- 
tem is the Willamowski-Rossler (WR) reaction-diffusion 
model [l2| given by 



ftc(r,t) = R[c(r,t)]+DV 2 c(r,t), 



(1) 



where Ri — K\C\ - 



R 2 = ^2Cic 2 — K-2C2 — n s c< 2 + and R3 = —n^c\Cz + 
ft_ 4 + ^503 — ^-503. These rate equations were derived 
from the mass action kinetics of a reaction scheme with 
quadratic kinetics where certain pool species are taken to 
be fixed. Here Cj(r, t) is the local concentration of species 
i at site r in a two-dimensional space of size L 2 with peri- 
odic boundary conditions. The parameters n±j are rate 
coefficients that contain the concentrations of the pool 
species that are fixed to maintain the system out of equi- 
librium. The diffusion coefficients D of all three species 
are taken to be equal. Even though the WR model is 
very simple, it exhibits a phenomenology with many 
features in common with those observed in chemical ex- 
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periments on the BZ reaction ^]J. Consequently, this 
model may be expected to capture the qualitative fea- 
tures of chemical systems whose chaotic attractors arise 
from a period-doubling cascade. The chaotic attractor 
for a certain set of rate coefficients is shown in the left 
panel of Fig. [I] 

A defect is characterized by its integer topological 
charge m top which is defined by Q 



j> V0(r, t) • dl = ±m t op, 



(2) 



where 0(r, t) is the local phase and the integral is taken 
along a closed curve surrounding the defect. A topolog- 
ical defect corresponds to a point in the medium where 
the local amplitude is zero and the phase is not defined. 
Typically only topological defects with m top = ±1 are 
observed. One-armed spiral waves with such defects at 
their centers are the only stable spiral waves for the 
CGLE 0. The phase has to be defined in order to 
apply the notion of defects to chaotic media. This is 
not a trivial issue and for many chaotic systems it is 
impossible to introduce a phase field. The rather sim- 
ple shape of the WR chaotic attractor (see left panel of 
Fig. 0} which arises from a period-doubling cascade ad- 
mits a simple definition of the phase and the WR model 
belongs to the class of chaotic-oscillatory media 0. We 
chose 0(r, t) = arctan[(c2(r, t) — C2)/(ci(r, t) — c?)] with 
( c i? c 2) = (8-0, 9.0) as the center of rotation. For such 
simple chaotic attractors, the particular choice of a phase 
variable for chaotic systems does not influence the results 
[l^ . In the center panel of Fig. d the phase field of the 
WR medium is shown for a certain set of parameters 
The topological defects can be identified as the termini 
of the white equiphase contour lines. As for the CGLE in 
the defect-mediated turbulent state 0, [2£j , the defects 
in the WR system for these parameter values rarely emit 
waves. They behave as passive objects and are merely 
advected by the surrounding chaotic fluctuations. The 
right panel in Fig.Hshows that the correlation function 
C(|Ar|) = (Re[A(r,t)A(r + Ar,t)]) rjt , where A(r,t) is 
the complex amplitude in the (ci, C2)-plane with respect 
to the center of rotation, decays exponentially with a very 
short characteristic length scale, verifying the existence 
of a turbulent state. This is the most typical character- 
istic for defect-mediated turbulence and has been found 
in the CGLE as well 

The fluctuations in the number of pairs of topolog- 
ical defects shown in the left panel of Fig. provide 
further evidence that a defect-mediated turbulent state 
can exist in chaotic-oscillatory media. The total num- 
ber of defects in the medium is exactly twice the num- 
ber of pairs because the net topological charge is con- 
served and equal to zero due to the periodic bound- 
ary conditions. Hence, topological defects can only be 
created and annihilated in pairs of opposite topologi- 
cal charge. One can easily derive a PDF p(n) for the 



number of defects provided that the defects are statisti- 
cally independent entities 0, • I n the stationary state 
and for periodic boundary conditions, the master equa- 
tion reduces to p(n) = p(n — l)c(n — l)/a(n) where c(n) 
and a(n) are the creation and annihilation rates, respec- 
tively. Provided that c(ri) = c = const and a(n) — an 2 , 
p(n) oc (c I a) n j (n!) 2 . The right panel of Fig. [2| shows 
that the PDF for the WR model agrees with the pre- 
dicted form of a squared Poissonian distribution reason- 
ably well. Moreover, the assumptions leading to this dis- 
tribution seem to be justified: the annihilation rate scales 
approximately with n 2 and the creation rate is approxi- 
mately independent of n as can be deduced from Fig. 03 

The fluctuations in the number of defects n do not 
have the properties one would expect if they arose 
from independent random events. If this were the 
case the power spectrum Sz,(/) = riniT^ool/2T | 

T 

J dt n(t) exp~* 27r ^ t | 2 for a system of size L would have 

— T 

a Lorentzian shape. The power spectrum does not have 
this form for the defect-mediated turbulence in the WR 
reaction-diffusion system. Figure El shows that Sl(/) oc 
l// 7 for intermediate frequencies with an exponent 7 
that is far from the value 7 = 2 expected for a Lorentzian 
shape. For k-2 = 0.072, we find 7 = 1.43 and, for 
k-2 = 0.075, 7 = 1.60. Although in both cases the sys- 
tem exhibits defect-mediated turbulence, the exponents 
are significantly different from each other. This implies 
that different chaotic-oscillatory media can have differ- 
ent second-order statistics. To confirm that these results 
are not specific to the WR system, we have also stud- 
ied the autocatalator reaction-diffusion system 0, [2^ 
which also has a chaotic attractor arising from a period- 
doubling cascade but with very different Lyapunov spec- 
tra. We find a power-law decay with values of 7 that 
are similar to those reported here for the WR model |23| . 
These results suggest that the power spectrum of n(t) 
may exhibit power-law decay with non-trivial exponents 
for defect-mediated turbulent states when the local dy- 
namics exhibits deterministic chaos. 

For large enough system sizes, oc L 2 for all /; 

thus, the power spectrum for large systems can be con- 
sidered to be the superposition of the power spectra of 
subsystems. Such behavior is expected in view of the 
short correlation length. This implies that the standard 
deviation E = ((n 2 ) — (n) 2 ) 1 / 2 scales as \J (n), which is 
proportional to L. Such a scaling is a consequence of 
the law of large numbers and has been observed for the 
CGLE as well. It follows that 7 and the low-frequency 
cutoff are independent of L. The low- frequency cutoff 
is due to the fact that n(t) is bounded. Its location de- 
pends on the density of defects in the medium. For lower 
densities the cutoff moves to lower frequencies. 

To understand how the non-trivial correlations in n{t) 
at intermediate time scales arise, we have analyzed the se- 
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FIG. 1: Left: Projection of the chaotic attractor in the (ci, C2)-plane of the homogeneous WR model for ki = 31.2, k~i — 
0.2, k,2 = 1-45, k-2 = 0.072, ^3 = 10.8, K-3 = 0.12, ka — 1.02, k-a = 0.01, ^5 = 16.5, K-5 = 0.5. Snapshots of the inhomogeneous 
c field in the defect turbulent state closely resemble this attractor. Center: Phase field for L = 128 Right: Correlation 

function in the defect- mediated turbulent state for L = 128 and different values of K-2. Note that n = 15.84 for K-2 = 0.072 
and ft = 8.81 for K- 2 = 0.075. 





15000 30000 5 10 15 20 25 
t n 



FIG. 2: Left: Time series of the number of defect pairs n for 
the parameters given in Fig. Qand L — 128. Right: Normal- 
ized histogram h(n) of n. The solid curve is the corresponding 
squared Poissonian distribution. 
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FIG. 3: Left: c(n) in the WR system for the parameters 
given in Fig. [I] and L = 128. The solid line corresponds to a 
constant creation rate. Right: a(n) in the WR system for the 
same parameters. The solid line corresponds to an increase 
in the annihilation rate proportional to n 2 . 
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FIG. 4: Power spectrum of the signal n(t) for different values 
of K-2 and L. For clarity, the curve for K-2 = 0.075 has been 
shifted down by one decade. The thick lines are to guide the 
eye. The thick solid line decays with 7 = 1.43 and the thick 
dotted line with 7 = 1.60. Note that the average period for 
one oscillation is T = 0.998 for «_ 2 = 0.072 and f = 0.980 
for k-2 = 0.075. 



ries of waiting times between consecutive creation events 
and consecutive annihilation events separately. In both 
cases, the waiting times are exponentially distributed and 
statistically uncorrelated as for a random walk. This im- 
plies that the correlations in n(t) are due to the interac- 
tion of creation and annihilation events. 

To compare our results for the power spectrum of n(t) 
in chaotic media with those in oscillatory media, we have 
simulated the CGLE [3,0 

d t A = A + (1 + ia)V 2 A - (1 + i(3)\A\ 2 A, (3) 

in a domain of size L 2 with periodic boundary condi- 
tions 24]. Here, A is the complex amplitude field and 
a and (3 are control parameters. The power spectrum 
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FIG. 5: Power spectrum of n(t) for different parameters in 
the CGLE and L — 128. The thick line is to guide the eye 
and decays with 7 = 1.9. Note that f = 12.7, 8.40, 5.43 from 
highest to lowest a. 



of the CGLE in the defect-mediated turbulent state has 
a power-law decay for intermediate frequencies and the 
scaling of Sl(/) with L, as well as the dependence of 
the low-frequency cutoff on the density of defects, agree 
with those in the WR medium. Yet, as Fig.Oshows, the 
exponent 7 = 1.9 which is very close to that expected 
for a Lorentzian form if the dynamics of n(t) is a simple 
random walk in a bounded domain. It is very different 
from the values observed for the chaotic media we ana- 
lyzed. Thus, different media can have different second- 
order statistics for n(t) depending on the underlying local 
dynamics. This argues against a universal description of 
defect-mediated turbulence. 

We have shown that defect-mediated turbulence can 
arise both in oscillatory and chaotic-oscillatory media 
and, thus, applies to a broad range of systems. While 
most common diagnostic measures do not allow one to 
distinguish between oscillatory and various chaotic me- 
dia, the fluctuations in the number of defects can be dif- 
ferent for different media: for the CGLE they resemble 
the form of a bounded random walk; for the WR and 
autocatalator models pronounced non-trivial correlations 
exist. Our results may be tested experimentally on sys- 
tems like the Belousov-Zhabotinsky reaction whose lo- 
cal temporal dynamics can be complex-periodic or even 
chaotic ^]|, an d where defect-mediated turbulence has 
been observed |(|. 
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